Macro- and micro-strain in GaN nanowires on Si(lll) 
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Abstract 

We analyze the strain state of GaN nanowire ensembles by x-ray diffraction. The nanowires 
are grown by molecular beam epitaxy on a Si(lll) substrate in a self-organized manner. On 
a macroscopic scale, the nanowires are found to be free of strain. However, coalescence of the 
nanowires results in micro-strain with a magnitude from ±(0.015)% to ±(0.03)%. This micro- 
strain contributes to the linewidth observed in low-temperature photoluminescence spectra. 
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I. INTRODUCTION 



Epitaxially grown semiconductor nanowires (NWs) have attracted large interest in re- 
cent years as means to fabricate nano-sized devices. They have been grown using various 
techniques like e.g. molecular beam epitaxy (MBE) and metal organic chemical vapor phase 
epitaxy. Often metal catalyst particles are used to initiate the growth. More recently, 
catalyst-free processes have been developed. The NWs promise to assist in overcoming 
many of the limitations of heteroepitaxy. Dislocations formed at the NW/substrate inter- 
face do not propagate along the N W axis, but remain at the interface or bend toward the 
free NW sidewalk. |2[ Residual strain is released elastically due to the extreme aspect ratio 



of typical semiconductor NWs. 3), |4j As a specific example, GaN NWs on Si are virtually 
free of extended defects despite the very large lattice and thermal mismatch between these 
materials. [5| Furthermore, the spectral position of the donor-bound exciton [(D^X^)] tran- 
sition in photoluminescence (PL) spectra of GaN NW ensembles is exactly the same as the 
one observed in strain free, bulk GaN. [6] In apparent contradiction, the linewidth of this 
transition exceeds values of 1 meV and is thus significantly above the value expected for ma- 
terial free of strain, jf^ 8] Recently, this minimum broadening invariably observed in the PL 
of GaN NW ensembles has been proposed to arise from the random position of donor sites 

nn 

and the corresponding energy distribution of excitons bound to these donors. |6|, [9] However, 
this intrinsic broadening mechanism cannot account for the fact that the linewidth observed 
for different samples may vary significantly, and may also approach values in excess of the 
energy separation of the donor-bound exciton transition in the bulk and close to the surface. 



Here, we investigate the morphological and structural properties of two GaN NW en- 
sembles on Si(lll) exhibiting a PL linewidth differing by more than a factor of two. Using 
x-ray diffraction (XRD), we determine the orientation distribution of the NW ensembles as 
well as their strain state. Position and shape of the x-ray diffraction peaks contain informa- 
tion about the strains and sizes as well as the orientation distribution of the NW ensemble. 
Both samples are found to be entirely relaxed on a macroscopic scale, but to possess a 
non-negligible micro-strain. 
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II. EXPERIMENT 



GaN NWs are grown on Si(lll) substrates by plasma-assisted MBE. Atomic N and Ga 
atoms are supplied by a plasma source operating at 500 W and an N 2 flux of 2 seem. Prior 
to NW growth, the Si substrate is exposed for 5 min to an atomic N flux for nitridation, 
which results in the formation of a 5 nm thick Si^N^ film. For sample 1 (2), GaN NWs were 
subsequently grown for 3 h at a substrate temperature of 780 (820V^C and N rich conditions, 
resulting in NW arrays with an average length of 1.0 (1.6) /xm.[lOJ High-resolution XRD 
measurements were performed at a temperature of 27°C using a Panalytical X-Pert PRO 
MRD™ system with a Ge(220) hybrid monochromator (CuKai radiation with a wavelength 
of A = 1.54056 A). The circles of the diffractometer are explained in Ref. The program 
Epitaxy™ was used for evaluation of the data. PL spectroscopy was performed using the 
325-nm line of a Kimmon He-Cd laser for excitation. The PL signal was dispersed by an 
80-cm Horiba Jobin-Yvon monochromator. The spectral resolution of the setup was about 
250 yueV. The investigations were done at a temperature of 10 K. TEM specimens were 
prepared by mechanical lapping and polishing, followed by Argon ion milling. TEM images 
were acquired with a JEOL 3010 microscope operating at 300 kV. 



III. RESULTS AND DISCUSSION 

Figures [T](a) and (b) display top- view scanning electron micrographs of samples 1 and 2, 
respectively. The azimuthal orientation of the NWs, as seen by the direction of their facets, 
is quite regular although the interfacial Si x N y layer is not expected to grow epitaxially. The 
NWs of sample 2 are distinctly thinner compared to those of sample 1. Coalescence of the 
NWs is obvious for both samples. 

Figure [2] shows near-bandgap PL spectra of the NW ensembles displayed in Fig. [TJ For 
both samples, the dominant line at 3.471 eV is due to the decay of the donor-bound exciton 



in strain-free GaN. 



• I2I, Q] Sample 1 exhibits a linewidth of 3.7 meV, more than twice larger 



than the one observed for sample 2 (1.6 meV). Broadening of bound exciton transitions 
by their energy dispersion as a function of the distance to the surface is expected to be 
more important for thin NWs, since the fraction of donors close to the surface increases 
quadratically with the diameter, j^, 0] In the present case, we would thus expect a narrower 



Figure 1. Scanning electron micrographs of (a) sample 1 and (b) sample 2 in top view. High- 
density arrays of facetted NWs are visible. The regular orientations of the facets point to a narrow 
distribution of the azimuthal orientation of the NWs. A significant fraction of the NWs constituting 
the ensembles have undergone coalescence. 

transition for sample 1 and not a broader one as observed experimentally. Consequently, the 
additional broadening observed for this sample has a different origin. As hypothesized in the 
following, a highly probable source of this broadening is microstrain caused by coalescence. 




Figure 2. (Color online) Near-bandgap PL spectra of samples 1 (thick line) and 2 (thin line) at 
10 K. The dominant line at 3.471 eV has a linewidth of 3.7 meV (1.6 meV) for sample 1 (2) and 
is attributed to the decay of the donor-bound exciton in strain-free GaN. 



4 



The two samples under investigation were fabricated at different substrate temperatures, 
which is known to affect both the nucleation density and the lateral growth rate. Lower 
temperatures as used for sample 1 result in a higher nucleation density as well as a more 
pronounced tendency towards lateral growth during nucleation, thus resulting in larger nu- 
clei. This scenario will inevitably lead to an enhanced coalescence in an early stage of NW 
growth, and significantly larger diameters of the final NWs as indeed observed in Fig. [TJ The 
misorientation of the NWs with respect to each other, as well as the dislocations accommo- 
dating this misorientation at the tilt or twist boundary, [ijj will introduce an inhomogeneous 
elastic distortion of the coalesced aggregate. Simultaneously, the strain introduced at the 
coalescence junction is more difficult to relax in the thicker NWs of sample 1. 4] All these 
arguments would point towards a larger contribution of inhomogeneous strain to the PL 
linewidth of sample 1, as also observed experimentally. In the following, we examine this 
hypothesis using high- resolution XRD. 

The out-of-plane orientation distribution of the NWs is determined from u scans near the 
GaN 00.4 reflection at several azimuthal orientations. The tilt of (2.3±0.2)° [(3.2±0.1)°] for 
sample 1 [2] are typical for GaN NWs grown directly on Si(lll). A pole figure of the GaN 
00.4 reflection (sample 1) shown in Fig. [3] displays a nearly isotropic angular distribution 
of the tilt angles of the GaN NWs analogous to a fibre texture. Figure H] demonstrates the 
dependence of the full widths at half maximum (FWHM) of the peaks of y?-scans (circles) 
performed using different reflections in skew geometry (sample 1) at different tilt angles \ 
(see 15|] for a sketch of the skew geometry) including in-plane scans for determination of 
the range of twist of the NWs. The full line shows the FWHMs extrapolated by: 



W hAl ( X ) = WV(WV360 + sin(x)) (1) 

where W^iix) * s the FWHM of a peak in the yj-scan for a reflection h, k, I with a given 
tilt angle x > an d Wqq is the FWHM of a peak in an in-plane scan revealing the true range 
of twist of the NWs. Often a <£>-scan can be measured in skew geometry more easily than in 
an in-plane geometry, then Wh,k,i(x) is measured and W 90 can be obtained from Eq. (Q]). 

The in-plane orientation distribution of the NWs as obtained from ip scans of the GaN 
10.0 reflection is about (2.3 ±0.3)° [(3.9 ±0.3)°] for sample 1 [2], and is again characteristic 
for GaN NWs on Si(lll). For both samples, coalescence boundaries will be the source of 
significant strain due to the broad orientation distribution of the NWs. This is illustrated 
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Figure 3. (Color online) X-ray pole figure of the GaN 00.4 reflection (sample 1) showing a nearly 
isotropic angular distribution of the tilt of the NWs. In this measurement x 1S restricted to 
0° < X < io°. 
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Figure 4. (Color online) Dependence of the full widths at half maximum (circles) of of the maxima 
of x-ray (/?-scans using different reflections (sample 1), each tilted by a certain angle x with respect 
to the sample surface. FWHMs extrapolated by Eq. ([!]) are shown as well (full line). 



in Figures [5] and [6j Fig. [5] shows a TEM micrograph of the interface region of sample 1, the 
GaN NWs on the Si x N y intermediate film. Strain contrasts are clearly distinguished near 
the bottom of the NWs and in regions where NWs meet and coalesce. Fig. [6] demonstrates a 
high-resolution TEM micrograph of the coalesced region of two GaN NWs (sample 1). The 
border acts as a small angle grain boundary and dislocations are formed in the vicinity to 
accommodate the tilt. The netplanes are locally curved due to the strain field of the defects. 



Figure 5. TEM micrograph of some of the GaN NWs on top of the Si x N y -film on Si 111 (sample 1). 
The local distribution of strain is visible from the contrasts on the micrograph. Those strain 
contrasts mainly due to bending of the NWs are found near the bottom of the NWs and near the 
coalesced regions. 
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Figure 6. High-resolution TEM micrograph of the coalesced region of two GaN NWs (sample 1). 
Dislocations are formed at the border similar to a small angle grain boundary and the netplanes 
are curved locally due to the strain field of the defects. 



Fig. [7] shows as another example a high-resolution TEM micrograph of an isolated crystal 
defect in a GaN NW (sample 1). This defect is probably caused by a stacking fault bounded 
by a dislocation. 

To determine the macroscopic strain state of the NWs, we first measured the lattice 
parameter of our Si substrates by triple-crystal diffractometry 111, [l6|. Using an analyzer 



crystal in front of the detector, the angular positions of the incident and the diffracted beams 



Figure 7. High-resolution TEM micrograph (Fourier filtered) of an isolated crystal defect (see black 
arrow) in a GaN NW (sample 1). This defect is probably caused by a stacking fault bounded by 
a dislocation. 



are determined precisely in order to obtain the Bragg angle 



'he necessary calculations for 



the substrate are performed in dynamical approximation [17|]. We thus obtained a value 
of 5.4299 A. Thanks to the small diameter of the NWs, the kinematical approximation is 
sufficient for the calculation of the angular position of their x-ray reflections. The positions 
of the symmetrical in-plane and out-of-plane reflections in u/2Q scans yield a = 3.1887 
(3.1885) A and c = 5.1855 (5.1854) A for the lattice parameters of sample 1 (2). Within 



16 



and 



18 



the error margin of ± 5 x 10~ 4 A, these values agree with those reported in Refs. 
for strain-free and pure GaN. This means, that our GaN NWs are indeed free of strain on 
average. 

The reciprocal space map shown in Fig. [8] for sample 1 confirms the fully relaxed state of 
the GaN NWs, because the diffraction vector h is pointing in radial direction from the GaN 
10.5 peak towards the origin of reciprocal space. The mosaic spread of the GaN NWs can 
be calculated from the width of the GaN peak perpendicular to the radial direction. The 
lateral correlation length is calculated from the reciprocal of the FWHM of the GaN 10.5 
peak measured parallel to the interface, i.e., along Q x . Here, rlu = \/2d is the reciprocal 
lattice unit, where d is the corresponding lattice plane distance. From the shape and the 
orientation of the GaN peak shown in Fig. [BJ we obtain a mosaic spread of (1.9 ± 0.2)° (i.e., 
close to the FWHM of the u scans) and a lateral correlation length of 76 nm (close to the 
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Figure 8. (Color online) Reciprocal space map near the GaN 10.5 and the Si 224 reflection for 
sample 1. The GaN lattice is fully relaxed, i.e., the reciprocal lattice vector of the diffracting planes 
h (arrow) is pointing to the origin of reciprocal space. rlu= X/(2d) stands for reciprocal lattice 
unit. 

actual lateral size of the NWs, see Fig. [1]). 

The widths of u/2Q scans across the GaN 00. n reflections with n= 2, 4, 6 increase linearly 
with the reflection order n, which is the result expected for broadening governed by micro- 
strain. The line shapes are nearly Lorentzian, so that we obtain the breadth (3f of the 
physically broadened profile j3f — (3h — (3 g where (3 g is the breadth of the apparatus function 
and (3h that of the measured curve [191 ] . Small asymmetries of the peaks are neglected here. 
We assume that each GaN NW is a coherently diffracting domain separated by air or a 
small-angle grain boundary in case of coalescence. These domains are tilted with respect to 
the interface and twisted around the surface normal by different angles reflected by broad 
peaks in the ui and ip scans. The u/2Q scans, instead, are broadened due to the finite size of 
the NWs and due to the inhomogeneous deformations inside the NWs. The integral breadth 
of the physically broadened line profile is /3/(2B) = A/I , where A is the integral intensity 
and Jo the peak intensity of the corresponding maximum. We assume that this breadth j3f 
of our sample consists of size and strain components (3 s and (3d- The size component (3s 
does not depend on the diffraction angle. The strain component (3d grows linearly with the 



order of reflection 



19]. We apply the reciprocal lattice notation: \d*\ = 1/d — 2sin(6)/A is 



the absolute value of the reciprocal lattice vector and (3% = /3/cos(9)/A. For a Lorentzian 
peakshape the two contributions have to be added [20] : 



P* f = P*+P D = (3* s + 2ed* 



(2) 
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Figure 9. (Color online) Williamson-Hall plot for the symmetrical GaN reflections of samples 1 (2), 
see dashed line (solid line). The inset shows the GaN 00.4 peak (cj/2@ scan) of sample 2 fitted 
with a Lorentzian curve. 



where e = fi* D /2d* = Ad/d is the strain. 20) We can interprete this strain as the variation of 
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d spacing within the NWs. The plot of (3j over d* is expected to give a straight line. 
Such a plot shown in Fig. is a Williamson-Hall plot in reciprocal-lattice representation. 
The intersection with the ordinate axis allows us to determine flg, corresponding to the size 
/ = l//3g. This is the size of an average domain in the direction of the diffraction vector, 
i.e., we essentially probe a length along the NWs. The characteristic size is I = 1.3 (0.3) fim 
for sample 1 (2). For sample 1 this length is comparable to the NW length, but for sample 2 
it is lower by a factor of five pointing to a smaller domain size inside the wires. The 
micro-strain e = ±2.8 x 10 -4 (±1.4 x 10~ 4 ) in the GaN NWs of sample 1 (2) is obtained 
from the slopes of the Williamson-Hall plots. These values describe the fluctuations of the 
lattice plane distances in the NWs, not their average. We estimate an error of ±1 x 10~ 4 
in the determination of the micro-strain. The inhomogeneous micro-strain is caused by 
residual defects in the GaN NWs like dislocations or even small angle boundaries arising 
after coalescence of the NWs. [l^ 

The important result obtained here is the fact that the micro-strain for sample 1 is 
indeed considerably larger than that for sample 2, just as the PL linewidths. A quantitative 
comparison of these values to the linewidth obtained in PL is, however, not straightforward. 
Coalescence presumably generates a rather complex strain field, including shear strain, while 
we only measure the ^-component of the complete strain tensor. To give an impression of 
the range of changes of the bandgap AE by a general strain with e zz = ±1 x 10~ 4 , we 
follow Ref. |22j and obtain AE^ = 1.65 meV and AE^ = 5 meV for pure biaxial (AEf,) and 
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pure hydrostatic (AEh) strain, respectively. Considering that at least part of the linewidth 
is caused by the random distribution of donors and the resulting energy dispersion, it is 
clear that only a fraction of the strain measured by XRD manifests itself in the linewidth 
observed in PL. We believe that this finding reflects the presence of dislocations in the 
coalescence junctions. Dislocations are known to act as nonradiative recombination centers 
and presumably suppress the PL from the junctions which act as the dominant source of 
strain measured by XRD. 

IV. SUMMARY 

We have performed structural characterization of MBE grown arrays of GaN NWs on a 
Si 111 substrate. The average lattice in our GaN NWs is fully relaxed. Typical residual 
strains in the wires can reach values from ±(0.015)% to ±(0.03)%. The local distribution of 
strain is visible in TEM micrographs. The observed micro-strain contributes to broadening 
of the lines measured by photoluminescence spectroscopy. 
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